import numpy as np
import pandas as pd
import statsmodels.api as sm
import math
import sys

### Operational audits and predicted race probabilitiess
dataBISG = pd.read_csv("/REDACTED/data/final/individualBISG2014_full_final.csv")
dataBISG= dataBISG[~dataBISG.predicted_prob_black.isna()]


##########################################################
#### TABLE A2
##########################################################
sys.stdout = open("/REDACTED/BIFSG_feature_count.txt", "w")

print(dataBISG.case.value_counts())
print(dataBISG.case.value_counts().sum())
print(dataBISG.case.value_counts(normalize = True))
print(dataBISG.case.value_counts(normalize = True).sum())
